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Introduction 


Forming a complement to our preceeding paper(J), we present some 
results for the first three harmonics derived from the simulation of 
diffusion-convection of galactic cosmic rays. We present also some de- 
pendences of the results to the modulation parameters. The results are 
discussed in comparison with observatipns and with the former studies of 
higher order anisotropy. 

In addition to this, we suggest the existence of the IMF- sense- 
dependent anisotropy of higher order which is discussed in detail in a 
separate paper(2). 


Theory and Model 


The cosmic ray anisotropy n(r, p) can be expressed as 
3 « 

r?(r ,p)M+£E ( (r , p )cos (m$' ) (r , p )sin (r«4> ' ) | PS(cos0' ) , (1) 

n-U =0 

where 0* and express the incident direction of cosmic rays, defined 


as( cf. Fig.l ) 

©'= 7C~0 , 7 1 + 

and 

rj , „ c (r,p)=(-l )"/ r S c (r,p)/F? c (r,p), 

i»;*(r,p)=(- 1 )"f'n s (r , p )/Fl c ( r , p ) . 

As shown in our preceeding paper (I), 
F$ c and F$ s .in Eq.(l) are given by 
the cosmic ray density U, stream S^, 
stress tensor T and heat flow 
tensor H ±j^ which are governed by 
the following equations. 



V- (CUV-k<”-VU)= -^(Ipy-VU), 
S=CUV~k (,) -vu, 

o 

Ti,-=-ic ( 2 ) i >/ 9 (VS )/ 9 . 

o 

II I )k = ~ K ^ i ikfgh ( VT ) f git . 


( 2 ) 

(3) 

(4) 

(5) 


Fig.l IMF-COORDINATE SYSTEM 

p ; Particle's momentum. 

V ; Solar wind velocity. 

®i; Unit vector in the direction 
of B in the 'away' sector. 
®3 = ei x V/ |ei x v| . 
e 2 = e 3 x e x . 
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where 

(VS ) ) --g (V 'S)6f gt (V -S)= -g^+S q Zqrr » 

and 

Zijk” £ x ; ’ € ^“ -7'kji • (6) 

Note that, in the IMF-coordinate system, Z±j^s in Eq.(6) represent the 
IMF-curvature and focusing. We first solve Eq.(2) in the model helio- 
sphere and, starting from the solution for U, we can calculate succes- 
sively Si, Tij and Hijjfc by Eqs.(3)-(5). In these calculations, the 
following scattering m.f.p. is assumed, 

(R/GV )exp } (l + aiCOS0fl (7) 

where Aq and ax are the parameters, R is the rigidity and is the polar 
angle in the heliocentric polar coordinate system. The calculations were 
carried out for various values of Aq and ax and, in this paper, the 
results in the following two cases will be presented. 

case I .... Aq - 0.016 a.u. and 0< = a\< = 3 9 

case II • • • • Aq = 0.032 a.u. and 0 < : a\<^3. 

It is noteworthy that, in Eqs.(3)-(5), Si, Tij, Hi jx having sub- 
script f 3 ! odd times are IMF- sens e-dependent and their associated anisot- 

ms 

ropies symbolized by r\ n with mark T s ! change their signs according to 
the alteration of the ’away* and * toward’ sectors whereas n n° with mark 
T c ? do not. In this paper, we restrict ourselves only to and the 

natures of mentioned above will be discussed in a separate paper (2), 

Results and Discussions 


The anisotropies n™ C?s in space can be expressed by the surface 
harmonics in the equatorial coordinate system and can be observed as 
solar daily variations at the Earth(3) . 

The resultant first space harmonic vectors for cases I and II are 
shown in Fig. 2, with black characters in the positive state and white 
characters in the negative state. It can be clearly seen that the 
harmonic vector in the positive state changes its phase toward earlier 
hours from that of 18h in the negative state, in accordance with the 
observed 22-year variation of the diurnal variation(4, 5, 6) . 

nh 


t SI ( to I l.m.c ) 

in=o 

R = 35.1 GV 


3 ] 



Fig. 2 FIRST SPACE HARMONIC VECTORS 

Arabic numerals represent the 
value of ax in Eq.(7). The black 
and white characters represent the 
positive and negative polarity 
states , respectively . 



Fig. 3 SECOND SPACE HARMONIC VECTORS 
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The resultant second space harmonic vectors are shown in Fig. 3. 

These vectors in both polarity states are mainly due to the n^ c -type 
anisotropy which is independent of 0 ( cf. Eq.(l) ). The magnitude 
pf) c of the anisotropy includes three terms expressing, respectively, 
the spatial derivatives of Si, the IMF- curvature effect and the IMF- 
focusing effect. On the basis of simulations, it is found that all these 
terms have significant contribution to the anisotropy and the IMF- 
focusing effect is not always most predominant. A similar independent 
anisotropy was obtained also by Bieber and Pomerantz(7) , on the basis of 
the following diffusion equation with respect to the pitch angle cosine 
( p=cos0 ) . 


a/oi.xi , t ) , ....a/pi ,x t , t ) 


at 


+ p V- 




'---2-1) 
dfi 


df (jitXj > t ) v 


dp 


2 L 


(1 


*) , 0 . (8) 


where Dyy is the Fokker-Planck coefficient for pitch angle scattering and 
L is the focusing length of IMF(S). They concluded that the anisotropy 
is principally a result of the focusing effect of IMF represented by the 
last term in Eq.(8). This anisotropy can be regarded as a special case 
which is applicable only for one dimensional diffusion along the IMF with 
infinite radius of curvature in the equatorial plane. In other words, 
this special anisotropy lacks the term expressing the IMF-curvature 
effect which has a dominant influence on the anisotropy as pointed out 
previously. 

Turning to the third order anisotropy, we obtained two dominant 
terms n§ c and r^ c . The eigen phases of the tri-diurnal variation arising 
from these two terms are almost orthogonal to each other, that is, 

(lh,5h) from n§ c and (3h,7h) from ni? It is noteworthy that one of the 
eigen phases (3h, 7h) arising from coincides with the observed phase 
around 7h in local time(9,70). Since the n^ c -type distribution is not 
symmetrical with respect to the field, this implies that the angular 



Fig. 4 THIRD SPACE HARMONIC VECTORS 



Fig. 5 RIGIDITY SPECTRA OF THE 
FIRST THREE HARMONICS 
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distribution of the third order anisotropy can not expressed generally 
in terms of only the pitch angle cosine y. In this respect, it is not 
suitable to discuss the third order anisotropy on the basis of the 
diffusion equation of Eq.(8). Fig. 4 shows the resultant tri-diurnal 
variations which are mainly composed of the above mentioned two types of 
distribution, i.e., qg 0 and n^ c . It is noted that the tri-diurnal 
variations obtained show polarity dependence. This dependence, however, 
is sensitive to the modulation parameters and, for instance, in case I 
we can get almost a polarity independent variation with increasing a^. 

Lastly, Fig. 5 shows the rigidity spectra of the first three harmon- 
ics obtained in case I for a^= 0 . The spectrum of the higher harmonic 
rises more steeply than that of the lower harmonic, and in the high 
rigidity region, the spectra of the second and third harmonic are almost 
proportional, respectively, to R and R 2 . This is due to X linearly 
increasing with rigidity R( cf. Eq.(7) ). Such a rigidity dependence is 
consistent with the observations ( 70, 7 7) . 
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